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Abstract 

A set of exactly computable orthonormal basis functions that are useful in 
computations involving constituent quarks is presented. These basis functions 
are distinguished by the property that they fall off algebraically in momen- 
tum space and can be exactly Fourier-Bessel transformed. The configuration 
space functions are associated Laguerre polynomials multiplied by an expo- 
nential weight, and their Fourier-Bessel transforms can be expressed in terms 
of Jacobi polynomials in A 2 /(k 2 + A 2 ). A simple model of a meson con- 
taining a confined quark-antiquark pair shows that this basis is much better 
at describing the high-momentum properties of the wave function than the 
harmonic-oscillator basis. 
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I. INTRODUCTION 



The harmonic oscillator basis has been used extensively in numerical calculations involv- 
ing confined constituent quarks 0. The advantage of this basis is that the basis functions are 
easily computed and Fourier- Bessel transformed. This permits the computation of matrix 
elements in either a coordinate or momentum representation. For problems involving light 
quarks, where confined quarks must be treated relativistically, matrix elements of kinetic 
energy operators that involve square roots, and matrix elements that involve momentum 
dependent Wigner and/or Melosh rotations are best computed in momentum space. Ma- 
trix elements of confining interactions, which are simple to compute in configuration space, 
require the evaluation of integrals over singular distributions fU in momentum space. The 
advantage of the harmonic oscillator basis is that it is straightforward to evaluate both types 
of matrix elements, completely avoiding difficult calculations. 

Although in principle it is possible perform very accurate model calculations using a 
harmonic oscillator basis, in practice the oscillator basis does not provide an efficient rep- 
resentation of meson eigenfunctions of mass operators that arise from a linear confinement 
plus one-gluon-exchange interaction. The reason for this is that the Coulomb part of the 
one-gluon-exchange interaction leads to momentum-space wave functions that fall off alge- 
braically rather than like Gaussians. This algebraic falloff is also consistent with predictions 
of asymptotic QCD. 

The purpose of this paper is to suggest the use of basis functions that have all of the 
advantages of the oscillator basis with the additional property that they fall off algebraically 
in momentum space. 

The configuration-space basis functions, sometimes called "Sturmians," 0, have been 
used successfully in atomic and chemical physics ||,|| calculations. In this paper, we exploit 
the fact that they have analytic Fourier-Bessel transforms, and these transforms fall off 
like polynomials in momentum space. This makes them well suited for models where it is 
advantageous to work in both configuration and momentum space. 
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The method is applied to low lying states of a constituent quark model of a meson. The 
model has a relativistic kinetic energy, a Coulomb interaction, a linear confining interaction, 
and a smeared out spin-spin interaction. Calculations of eigenfunctions and eigenvalues 
of this constituent quark mass operator, obtained by projecting the mass operator on the 
subspace of the Hilbert space spanned by a finite number of these basis states, give converged 
eigenfunctions and good variational bounds on the mass eigenvalues. In particular, the high- 
momentum tail of the wave function is described quite well in a truncated basis. 

We also discuss briefly the extension of the method to cases where the configuration-space 
wave function has an integrable fractional power singularity at the origin f|[7|]. 



In configuration space the radial basis functions are given in terms of polynomials in r 
times an exponential. An orthonormal basis set is then given by || 



II. ELEMENTARY CONSIDERATIONS 
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where x = Ar, L^ l+2 \x) is the associated Laguerre polynomial 




(2) 



with a = 21 + 2 and 



N ln = A- 3 (i) 



T(n + a + 1) 



(3) 



is defined so that the (f> n i(r) satisfy the orthonormality condition 




(4) 



The momentum space functions are the Fourier-Bessel transforms of the n z(r)'s: 




(5) 



They satisfy the orthonormality condition 
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S nm = / k'dk^k)^)- (6) 
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With these definitions both <f> n i(r) and (f> n i(k) are real. 

These functions, along with spherical harmonics, can be used to expand either the mo- 
mentum or coordinate representation of approximate eigenstates. The coefficients in the 
expansions 

) = c mm<Pni{r)Yi m (r) (7) 



nlm 



and 



are related by 



^(k) = C nlrn<fird(k)Yi m (k) (8) 
nlm 



Cnlm ( ^) C n i m . (9) 



These relations follow directly from the formula for the spherical expansions of plane waves 
||. The phase factor ensures that ip*(k) = k) for real "0( r )- 

We now consider the momentum-space wave function, given by the Fourier-Bessel trans- 
form of ([I]): 

Mk) = \ - / r 2 3l (kr)Mr)dr (10) 

V 7T JO 

First, we note that the Fourier-Bessel transform of the unnormalized function r l e~ Ar is 

Second, we note that multiplication of the unnormalized function r l e~ Ar by r is equivalent 
to the operation —d/dA on its Fourier-Bessel transform: 

^ r J,,, - m„ -a,_ 9 /'72(2A)(2* ; ) i (; + l)!' 



Wo ''^''> r+e --aA^^AW^J' (12) 

Applying —d/dA to the right hand side of Eq. flT2"D n-times gives the original function 
multiplied by a polynomial of degree n in r = l/(/c 2 + A 2 ). In addition the coefficient of 



the leading power of r is of the same sign as the coefficient of the leading power of r in 
the associated Laguerre polynomial. The Fourier Bessel transform, being unitary, preserves 
orthonormality. Thus the Fourier-Bessel transform of <p n i has the form 

~ W = (A2 + % +2 >< Qn(r), (13) 



■>nl\ 



where the Q n (r), are orthogonal polynomials in r with weight k 2l+2 / (A 2 +k 2 ) 2l+4 . To identify 
Eq. (|13D with the Fourier Bessel transformation of Eq. ([[J) it is enough to choose the phase 
of the normalization constant so the coefficient of r n in Q n {r) has the same sign as the 
coefficient of r n in L^ +Z (r). 
Specifically, one can write 

Mk) = ^ {y2 y l iy+2 K n (u), (14) 



'N ln 

where y = k/A and K n (u) is a polynomial in u = l/{y 2 + 1) which satisfies the following 
orthonormality condition: 
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8nm = / y 2 dy 
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K n {u)K m {u). (15) 



_{y 2 + l) l + 2 

The integral in Eq. fll5|) can be transformed by the following variable change [ IU[ : 



il + x . dx l-x 

y = \lTzr^> dy = ~ — —> u = ^r~ ' ( 16 ) 



(1 - X) 2(1 + X) 2 

to the integral 

5 nm = 2~ 21 - 4 J 1 dx(l - x) l+ l(l + x) l+l 2K n (u)K m (u). (17) 
The polynomials K n (u) can be expressed in terms of Jacobi polynomials: 

K n (^)=Pn +2 ' l+l2 \x), (18) 



where 



pafi( T \- T (* + n + l) A (n\ T(a + + n + m+l) 
n {X) n \T(a + P + n + l)^ [m) 2™T(a + m + 1) 1 > ' { J) 



This gives the following expression for <p n i{k): 



where the normalization constant N n i is 

a 3 r(n + / + |)r(n + / + |) 

' n 2(2n + 2Z + 3) n!r(n + 2Z + 3) ' 1 ' 

At this point the normalization is determined up to an overall phase. The phase is fixed 
by requiring the sign of the leading power of r in the associated Laguerre polynomial be the 
same as the sign of the leading power of 1/(A 2 + k 2 ) appearing in the Jacobi polynomial. 
Inspection of equations @, fll9f) , and ( ^0|) show that the phases are (— ) n in both cases. This 
shows that ( p0|) is the Fourier Bessel transform of (p. 



III. APPLICATION 

To test the method, we consider a model |TT| of a meson consisting of a bound quark- 
antiquark pair. The mass operator (Hamiltonian in the center-of-momentum frame) is 

M = 2Vk2 + w? - ^ + (3r + 7 + 2s ( s + X ] ~ 3 (22 ) 

The parameter values a s = 0.5, /3 = 0.197 GeV 2 , 7 = 0.77 GeV, quark mass m = 
0.36 GeV, and tq = .66 GeV^ 1 . The parameters 7 and p are chosen to provide a fit to the 
physical 7r and p masses of 0.140 GeV and 0.784 GeV, respectively. 

The mass operator is diagonalized using 10, 20, 40 and 80 basis states. The momentum 
scale A = 2.0 GeV is chosen to minimize the mass eigenvalue for 10 basis states. The 
eigenvalues are given in Table |. The corresponding momentum wave functions are shown 
in Fig. 0. Using only 10 basis states, the wave function is quite stable up to k = 15 GeV. 

The calculations are done using recursion relations to compute both the associated La- 
guerre and Jacobi polynomials: 

in + l)4+i(>) = (2n + a + 1 - x)L^(x) - {n + ajL^x), (23) 
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and 

2(n + l)(n + a + (3 + l)(2n + a + 0)p£f (24) 
= [(2n + a + /3 + l)(a 2 -/3 2 ) 

+ (2n + a + /5)(2n + a + (3 + l)(2n + a + /3 + 2)a;]Pi Q ' /3) (x) 
-2(n + a) (ti + /3)(2n + a + /3 + 2)P^f } (x). 

For comparison purposes, we have also solved the eigenvalue problem with a harmonic 
oscillator basis, with 

^k) = ^=y l ti +1 *\y 2 )e-y 2 / 2 . (25) 

Vn 

The mass operator is again diagonalized using 10, 20, 40 and 80 basis states. The momentum 
scale A = 0.9 GeV is chosen to minimize the mass eigenvalue for 10 basis states. The 
eigenvalues are given in Table [n[ Convergence to the exact eigenvalue is much slower than 
with the polynomial basis. The corresponding momentum wave functions are shown in 
Fig. |2|. For each basis size, the momentum wave function is stable until a critical momentum 
is reached, after which it exhibits a Gaussian falloff characteristic of a truncated oscillator 
basis. For 10 basis states, this cutoff is about 5 GeV, which may be adequate for use in 
calculations involving moderate momentum scales (say, 1 GeV or less), but much better wave 
functions can be obtained with the same computational investment using the polynomial 
basis. 



IV. INTEGRABLE SINGULARITY 



The exact / = configuration-space wave function for the model mass operator of Eq. (p0| 
has an integrable singularity [[|[7| : 



^{r) ~ ± (26) 



as r — > 0. This singularity is seen in numerical calculations based on B-splines in \Y1\ and 
persists in the presence of the confining interaction. 



Naive application of this wave function can lead to unphysical predictions in calculations, 
such as decay widths, which depend on the value of the wave function at the origin. This is 
a limitation of this naive model which must be cured by "correcting" the form of the mass 
operator. Independent of these considerations, it is possible to generalize the methods of 
this paper to include this mild singularity in the basis functions. 
The relevant configuration space basis functions (for I = 0) are 

1 e~ Ar 



<f>n(r) 



/AC r a 



-L 2 n {1 - a) (2Ar 



with 



T(3-2a + n) 



(27) 



(28) 



n! (2A)3-2a 

The Fourier Bessel transforms can be computed explicitly by observing that the Fourier 
Bessel Transform of 



£ n (r) := r 



(29) 



is 



2 1 r (n + 2-a) . .. n , -i,A\i 
~TTTn , q \ / ^ — r^- sin \[n + 2 — a) tan (-) . 



which follows from the expressions ||l3f 

J e- Kr r»J v {kr) = T(n + v+ l)r~ i,x+1) P~ u (A/VP + A 2 



and 114 



^ (MW)= J j#im , 

V * (/< + 1/2),/5^ 



It follows that 



<f> n (k) 



(30) 



(31) 



(32) 



I 2 n 



^Nn V vr m=0 



f n + 2 - 2cA (2A) m T(m + 2 
v n-m J m! (A 2 + k 2 f m+2 - a )l 2 



a) sin[(m + 2-a)tan- 1 (|)] ^ 



The relevant basis functions and their Fourier Bessel transforms can also be computed 
for / using the same methods, but these solutions are not singular at the origin for I > 0. 
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V. CONCLUSION 



We have presented a set of exactly computable orthonormal basis functions that are useful 
in computation involving constituent quarks. The basis functions can be computed easily in 
both position and momentum space, making it simple to calculate matrix elements in either 
space. A simple meson model shows that this basis is much better at describing the high- 
momentum properties of hadronic wave functions than the usual harmonic-oscillator basis, 
even when only a few configurations are used. The utility of this basis is not diminished in 
the presence of a spin-spin interaction of a scale to produce the physical n-p mass splitting. 
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TABLES 



TABLE I. Eigenvalues of the pion mass operator using the polynomial basis. 



basis 


M 


10 


0.14071 


20 


0.14037 


40 


0.14019 


80 


0.14013 


TABLE II. Eigenvalues of the pion mass operator using the oscillator basis. 


basis 


M 


10 


0.14826 


20 


0.14361 


40 


0.14186 


80 


0.14109 
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FIGURES 

FIG. 1. Pion momentum wave function using the polynomial basis. The dotted, dot-dashed, 
dashed and solid curves correspond to using 10, 20, 40 and 80 basis functions, respectively. 

FIG. 2. Pion momentum wave function using the oscillator basis. The dotted, dot-dashed, 
dashed and solid curves correspond to using 10, 20, 40 and 80 basis functions, respectively. 
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